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Abstract 

We have calculated the binding energy of the hybrid quarkonium up to the order 
Cn _ 5 

of C(l/mg) and found their decay constants scale like rriq as mg — * oo. The 
and ++ hybrid quarkonium is exactly degenerate in the limit txiq — ► oo while the 
■^J- ■ 0(\/mo) correction renders the mass lower than that for ++ . The 1 ^ and 

1 H hybrid is nearly degenerate and lies 0.7 GeV lower than the . If we use 
nib = 4.8GeV the mass of the 1 ^ bbg state is (10.75 db 0.20) GeV, which is the 
ON ' lightest exotic hybrid quarkonium. 
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Quantum Chromodynamics (QCD) is believed to be the correct theory of strong in- 
teraction. In the massless light quark limit, QCD has exact SUl(3) x SUr(3) symmetry. 
On the other hand, QCD has spin-flavor symmetry in the heavy quark limit tjiq — > oo. 
In both cases we have an effective theory: the chiral perturbation theory (%PT) and the 
heavy quark effective theory (HQET). xPT is very useful for the low energy dynamics of 
strong interaction, while HQET provides a consistent framework to study hadrons con- 
taining one heavy quark. Symmetry plays a fundamental role in particle physics. It is 
important to explore novel features of the exotic system (QQg) with two heavy quarks in 
the limit niq — ► oo. 

From quark model we know that a qq meson with orbital angular momentum 1 and 
total spin s must have P = (— 1)' +1 and C = (— l) l+s . Thus a resonance with J PC = 
, + ~, 1~ + , 2 + ~, • • • must be exotic. Such a state could be a multiquark state or gluonic 
excitation such as hybrids, glueballs. Very recently there appears experimental evidence 
for a light J PC = 1~ + exotic 0, ^ |||, f|, || . The emergence of evidence for hybrids indicates 
the presence of dynamical glue in QCD. 

The traditional theoretical approach to the gluonic excitations in mesons falls into 
two classes: the string or flux tube model and the constituent glue model with the glue 
confined by a bag or potential. In contrast, QCD sum rule (QSR) starts from the first 
principle like the lattice gauge theory and incorporates the nonperturbative effects via 
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various condensates in QCD vacuum. It has proven successful in extracting low-lying 
hadron masses and decay constants |J. Several collaborations [0-0 have studied the 
masses of light hybrids with QCD sum rules. The mass of heavy hybrid quarkonium has 
been studied by one group PH|. But unfortunately the QCD sum rules for the exotic {QQg) 



states with J PC = , 1 + are not stable, which prevents a reliable extraction of their 
|T0| . Recently the masses and decay widths of hybrid mesons containing one heavy 



masses 



quark have been studied in the framework of HQET |TTJ . 

In this work we shall investigate the binding energy of hybrid quarkonium (QQg), 
the scaling behavior of their decay constants in the limit of tuq — > oo and their I/tuq 
corrections. As we shall show later, all of our sum rules are stable after the heavy quark 
mass is separated out. 

The interpolating current for the J PC = l~ + ,0 ++ heavy hybrid quarkonium reads 



J„(x) = Q(x)g s YG%(x)^Q(x) , (1) 



and for the J PC = 1+" 0" 



ones 



Jj(x) = q{x)g s ^^Gl u {x)^h v {x) . (2) 



Denote the J PC = 0++, 0~ hybrid mesons by Hi and the decay constants by 
fi, i = 1,2, 3, 4. For example, fi is defined as: 

(01^(0)1^)=/^, (3) 

where is the H\ polarization vector. 
We consider the correlators 

i f Ae^(o|T{j,W, 4(o)}|o> = -{g,u - ^)Mp 2 ) + ^Mp 2 ) , (4) 

J p p 

1 1 d 4 xe* kx (0\T{Jl(x), J 5 J(0)}\0) = -far - ^)n 3 (p 2 ) + ^n 4 ( P 2 ) . (5) 

The imaginary parts of IIx, n 2 , n 3 , n 4 receive contributions from the Hi, H 2 , H 3 , H± hybrid 
intermediate states respectively. 

The dispersion relation for Ilj(p 2 ) reads 



rV) = / p{s J . ds , (6) 

J s — p l — le 

where p(s) is the spectral density. 
At the phenomenological side 

n(p 2 ) = T7o q + excited states + continuum . (7) 

M 2 — p l 

In order to derive the spectral density in full QCD we can either use the Cutkosky cut- 
ting rule or employ the dispersion relation and previous results of QCD sum rule for heavy 
quarkonium ||12j| . With both methods we have obtained the same results independently. 



For example, the spectral density in full QCD for H\ reads 



Pi(s) = —^ : \1 ^(1+ *y J -i '-dx+ ^— J \/l ^(1+ 5t 



(8) 
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which agrees with the expressions in Jl0| 



Note the gluon condensate and the perturbative term is of the same order in the limit 
of rriQ — > oo, which enables us to explore the scaling behavior. Let M = 2mg + A, where A 
is the leading order binding energy. Expanding (H), (M), (g) to the leading order of 1/rnq 
and making Borel transformation to suppress the continuum contribution, we arrive at 

-f 2 e t = / m^p (e)e ?de , (9) 



4m Q jo 

where p (e) = < 0|# 2 G 2 |0 > Je, a s = 2 B /2 - and < 0|^G 2 |0 >= 

3 ' QCD 

0.48 GeV 4 . The coefficients cij etc are listed in TABLE I. E c is the continuum threshold. 
Starting from E c we have modeled the phenomenological spectral density with the parton- 
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like one mqp^e). The triple gluon condensate is suppressed by a factor of l/m^ so it 
does not contribute up to the order 0{l/mQ). From @ we can introduce a scaled decay 

_ 5 

constant F = rag* f, which remain constant in the limit tuq — > oo. 

We further expand the hybrid mass M and decay constant / to the order of 1/rriQ, 

M = 2m Q + A + / = m Q (F + ^), A 1} F 1 is the 0(l/m Q ) correction to A and F 
respectively. Similarly we have, 

i F 2 ( \ -4- —\ i-E 

L\2F X F 1 l J 4 = / Pi(e)e-^de , (10) 



4 L T Jo 

where Pl (e) = Qffe 5 v^+ £ < 0|^ 2 G 2 |0 > e^. 
With (g) and (0) we get 



a _ Jo c spo(s)e ?ds 
Jo Po(s)e rrfs 



F 2 A 2 r E <= 

— (Ai + —)= Pi(e)(e - A)e'^de , (12) 



-A 



4 v 4 7 7o 

1 p l {e){l + — r -)e-—de. (13) 



2 Jo r v /v T 
Requiring the absolute value of the gluon condensate be less than 30% of the leading 
perturbative term with continuum subtracted we get the lower limit of the Borel param- 
eter T and the continuum threshold E c . Typically in our analysis the gluon condensate 
contribution in the whole sum rule is less than 25% for 1 h and 10% for the other three 
states starting from T > 0.9 GeV. We have kept four active flavors and let Aqqj^ = 220 
MeV. Varying Aq CD from 220 MeV to 300 MeV the final result changes within 5%. We 
want to emphasize that all the sum rules for A, A\,F, F\ are stable in the large interval 
0.9 < T < 3.0 GeV. 

As a consistency check we have also fitted the left and right hand side of Eq. @ 
directly with the most suitable parameters A, F, E c directly in the working region of the 
Borel parameter. The results agree very well with those derived from the derivative method. 

We collect our final results in TABLE II. For A there is an error about (±0.2) GeV. 
For the exotic state Hi we present the variations of A with T and E c in FIG. 1. We also 
show the left and right hand side of (0) with the central values of A, F, E c in TABLE II for 
Hi. Both sides agree within one percent in the region 0.5 < T < 3.0 GeV as can be seen 
from FIG. 2. 
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In the limit irtq — ► oo the and ++ hybrid quarkonium is exactly degenerate while 
the O{l/mo) correction makes mass lower than that for ++ . The 1 h and l + ~ 
hybrid is nearly degenerate and the 0(1/ mo) correction does not split their masses. The 
1~ + hybrid lies 0.7 GeV lower than the hybrid. The mass of the lowest exotic hybrid 
quarkonium bbg is (10.75 ± 0.20) GeV if we use m b = 4.8GeV. 

In summary we have calculated the binding energy of the hybrid quarkonium up to the 
order of 0(1/ mo). We have found that the decay constants of the hybrid mesons scale 
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like ttlq as niq — > oo. Moreover, the scale of the binding energy is solely set by the gluon 
condensate. In the sum rules for the light hybrids the four quark condensate and three 
gluon condensate are not precisely known, which renders the extraction of the light hybrid 
mass rather difficult [^, |[. In our calculation of hybrid quarkonium masses we have not 
considered the error due to the bottom quark mass. Another possible source of error is the 
lack of E c stability. 
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Figure Captions 



FIG. 1 The variations of A with T and E c for H\. T is in unit of GeV. 

FIG. 2 The variation of the right and left hand side of Eq. @ with T is plotted 
as solid and dotted curves respectively for H± with the values in TABLE II. 



TABLE I. The coefficients etc in po( e ) an d 





1-+ 


0++ 


1+- 


0~ 


a 
b 


256 


384 


128 


256 


945 
1 

24 
128 


945 
1 

16 

3648 


405 
1 

144 
1216 


945 
1 

24 
2176 


c 
d 


945 
1 


10395 


10395 
7 


10395 


144 


_ 96 


288 


_ 48 



TABLE II. The values of A,Ai,F,Fi,E c for the hybrid quarkonium. A, E c is unit of 
GeV, Ai is unit of GeV 2 , F is in unit of GeV~ and F\ is in unit of GeV~. 





1-+ 


0++ 


1+- 


o— ■ 


A 


1.22 


1.97 


1.21 


1.97 


Ai 


-0.44 


-0.44 


-0.42 


-0.92 


F 


0.244 


0.89 


0.208 


0.515 


Fx 


-0.064 


-0.105 


-0.036 


-0.42 


E c 


1.6 


2.3 


1.5 


2.3 
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